Spin splitting in graphene studied by means of tilted magnetic-field experiments 
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We have measured the spin splitting in single-layer and bilayer graphene by means of tilted 
magnetic field experiments. Applying the Lifshitz-Kosevich formula for the spin-induced decrease 
of the Shubnikov de Haas amplitudes with increasing tilt angle we directly determine the product 
between the carrier cyclotron mass m* and the effective g-factor g* as a function of the charge 
carrier concentration. Using the cyclotron mass for a single-layer and a bilayer graphene we find an 
enhanced (/-factor g* — 2.7 ± 0.2 for both systems. 
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The half-integer quantum Hall effect in single-layer 
graphene (SLG) [TJ [5] and the unconventional quantum 
Hall effect in bilayer graphene (BLG) |3j reveal spin- and 
valley-degenerate relativistic Landau levels. Due to the 
extremely large Landau-level splitting [H [5] , completely 
resolved levels can be observed up to room temperature 
[BJ. However, even at very high perpendicular magnetic 
fields the Zeeman splitting within one Landau-level is 
negligible smaller compared to the Landau-level splitting 
and, more importantly, the Landau-level width generally 
exceeds the spin-splitting. Exceptionally, the zeroth Lan- 
dau level in SLG becomes extremely narrow at magnetic 
fields B > 20 T [3], which allows an experimental obser- 
vation of a spin-related gap opening at magnetic fields 
B > 20 T 7\ . Another observation of a spin degeneracy 
lifting with an effective ^-factor g* — 2 was reported for 
v = ±4, in SLG for magnetic fields B > 30 T, combined 
with lifting the valley-degeneracy at v = ±1 [8 . 

In this paper we determine the spin splitting of broad- 
ened Landau levels for SLG and BLG by measuring 
Shubnikov-de Haas (SdH) oscillations in tilted magnetic 
fields. This technique allows adjusting the ratio between 
the spin splitting and the Landau level splitting, by con- 
trolling the ratio between a total magnetic field and a 
component perpendicular to a two-dimensional graphene 
flake. Using the well-established Lifshitz-Kosevich for- 
mula [9, 10 we determine the product of effective g-factor 
and cyclotron mass, m*g* , from the angular dependence 
of the SdH amplitudes and we find that g* is enhanced 
compared to the free electron value. 

We have fabricated field-effect transistors from SLG 
and BLG, by micromechanically exfoliating graphene 
flakes from graphite. The flakes were deposited on top of 
a Si/Si02 wafer, structured into a Hall-bar and covered 
with Au/Ti contacts [11] . Charge carriers are introduced 
by applying a gate voltage on the conducting Si substrate. 

We present a detailed analysis on the spin splitting in 



a SLG sample made from Kish graphite with a mobil- 
ity /i = 0.8 Vm _2 s _1 and BLG sample originating from 
natural graphite with a mobility /i = 0.3 Vm _2 s _1 . Two 
other devices, one SLG and one BLG sample, showed 
qualitatively similar results. 

To determine the spin-splitting we have measured the 
longitudinal resistances R xx as a function of charge car- 
rier concentration n at a constant perpendicular mag- 
netic field. We adjusted the total magnetic field B tot for 
each tilt angle such that the normal component B n is 
the same (inset to FigfTI). The value of B n was verified 
by measuring the Hall resistance of the devices in the 
non-quantized regime. 

In Fig. [I] we show the experimental R xx (n) dependen- 
cies for SLG at B n = 6 T (a) and for BLG at B n = 8 T 
(b). R xx shows Shubnikov-de Haas oscillations with max- 
ima whenever the Fermi energy is situated in the mid- 
dle of a spin- and valley-degenerated Landau level En, 
N = 0,1,2,... being the Landau-level index. For the 
higher Landau levels (N > 2) the longitudinal resis- 
tances do not exhibit zero minima indicating that the 
level broadening is comparable to the cyclotron energy 
at these perpendicular magnetic fields. 

When increasing B to t at a constant B n the oscillation 
amplitudes for both BLG and SLG are reduced. From 
this reduction we determined the spin-splitting. We use 
the Lifshitz-Kosevich formula for systems with a gen- 
eral dispersion and we specifically include spin-splitting 
[TO] with an effective g-factor g* QH [T3] and tilted 
magnetic fields [T3] . The oscillatory contribution to the 
longitudinal resistance can be described as [2]: 
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where S(E)\ E=E is a extremal cross section of the Lan- 
dau orbits in the fc-space, A is the oscillation ampli- 
tude and ipB is Berry phase, tp B = 7r for SLG [TJ |2], 
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FIG. 1: (color online) Shubnikov de Haas oscillations in SLG 
(a) at T = 1.3 K and in BLG(b) at T = 0.4 K as a function 
of the carrier concentration for different total fields B to t or 
tilt angles 6, respectively. When varying 9 the total field 
Btot is adjusted such that the perpendicular field B n remains 
constant, i.e. Btot = B n / cos 8. The oscillation maxima are 
marked with the corresponding Landau level numbers TV. The 
inset schematically shows this tilting configuration. 



ifB = 27r for BLG [3 . The amplitude A contains a 
monotonic n-dependent part, a temperature dependence, 
a 73 n -dependent contribution and a damping factor due 
to spin splitting depending on the total field B tot . At a 
constant temperature and perpendicular magnetic field 
this -B to t-dependence of the SdH amplitude A for charge 
carriers with cyclotron mass m* and effective ^-factor g* 
is given bv [T21IT4] : 



A = A (iV)cos 
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with cyclotron mass pQ: 
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and A (N) is constant for a given N. 

For the spherical Fermi surface in SLG and BLG with 
a Fermi wave-vector Icf = y/nn, the extremal cross sec- 
tion of the Landau orbits is S(E)\ E=E = irkp = mr 2 
and Eq. |l]) yields the concentration-dependent resistance 
oscillations as we observe them in our experiments: 

R xx = ^Icos f —- n + n + Lp B J = A cos (-v + n + Lp B ) , 

(4 ) 




FIG. 2: (color online) Normalized oscillation amplitudes as a 
function of total field Btot at a constant perpendicular field 
B n in SLG (a) and BLG (b). Error bars represent standard 
least squares fitting errors in the determination of A. Solid 
lines are fits to Eq. [2] with m* g* as a fit parameter. 



where v = (hn)/(eB n ) is the filling factor. As expected, 
the oscillation period, (2eB n )/(hn), is independent on 
the band structure of the 2D material and only depends 
on the filling factor. 

To accurately determine the experimental oscillation 
amplitudes we have fitted our experimental data R xx (n) 
to Eq. [4] in two steps. First we determined the oscillation 
period and a smooth background using all oscillations 
measured for a wide range of the carrier concentrations. 
Second we fitted the oscillation amplitudes A for each 
individual oscillation using the above determined period 
and background as fixed parameters. In Fig. [5] we show 
the final results of this fitting procedure for the SdH am- 
plitude as a function of the total magnetic field for dif- 
ferent Landau levels N. For clarity all amplitudes are 
normalized to Aq. 

The experimentally observed reduction of the SdH am- 
plitudes can be qualitatively visualized in a simple den- 
sity of states (DOS) picture of a Landau level as depicted 
in Fig. [3^,. In a purely perpendicular magnetic field the 
Landau level width exceeds the spin splitting and the 
DOS of the spin-down state (orange, horizontally dashed 
in Fig. |3k) overlaps with the one of the spin-up states 
(red, vertically dashed) to one broad Landau level. When 
increasing B tot by leaving B n constant, these two states 
move apart yielding an additional broadening of the Lan- 
dau level with a reduced DOS in the center (green, solid 
areas in Fig. pa). Eventually, when the spin splitting ex- 
ceeds the level width a minimum between two distinct 
levels starts to develop in the DOS. This scenario is in- 
deed observed experimentally in SLG (Fig. [3b). The SdH 
maxima corresponding to the N = 9 and N = 10 Lan- 
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FIG. 3: (color online) Schematic representation of the density 
of states for a Landau level with an increasing total magnetic 
field Btot (from the bottom to the top) at a constant perpen- 
dicular component B n (a). Panel (b) shows this scenario as 
measured experimentally for the N = 9, 10 maximum in SLG 
at a constant perpendicular magnetic field B n =5T. 



dau levels at B tot = B n =5 T do not show any splitting. 
Increasing of the total field at a constant perpendicular 
component leads to a reduction of the oscillation ampli- 
tude and eventually appearance of spin-resolved peaks 
at the highest field of 28 T. However, this splitting is 
not yet enough to determine the energy difference by e.g. 
activation measurements. 

A quantitative analysis of this decrease of the SdH am- 
plitudes with increasing total magnetic field is done by 
fitting the data to Eq. (pi) with m*g* as a fitting param- 
eter (solid lines in Fig.[2F. The values for m*g* obtained 
are plotted as a function of the charge carrier concentra- 
tion in Fig.[4]for SLG (a) and BLG (b). 

For both SLG and BLG the product m*g* increases 
with concentration, which can be mainly attributed to 
the concentration dependent cyclotron mass m* of par- 
ticles with a linear pQ and hyperbolic dispersion [15] as 
predicted by Eq. [3j 

The dashed lines in Fig. |4^ show the calculated de- 
pendence of m* g* for g* — 2 and g* = 2.7 using 
m*{n) = (h/c) yplra [lj. The shadowed areas represent 
a 10% uncertainty of this calculation mainly due to the 
experimental errors and some uncertainty in the Fermi 
velocity [16]. 

For SLG (Fig.[4k), the increase of m*g* with n is sym- 
metric for electrons and holes (i.e. negative and positive 
n in the figure). A best fit using m*(n) for SLG yields 
g* = 2.7±0.2 (the error is the standard deviation). This 
finding is shown directly in the inset of Fig. [4k, where we 
plot the value of g* determined in the middle of each Lan- 
dau level N for different perpendicular fields B n . Within 
an experimental error g* does not show any dependence 
on N or B n . 

For BLG (Fig. Eb) the experimental situation is more 
complex as the observed increase of m*g* with n is 



E 



0.18 


-4 


-2 2 


4 








12 


4 
2 




— 9*=2.7- 

\ ® 


0.06 
n nn 


- /I "^^\ M,''o B=5T - 

g =2 ^ \ / ^ 

Inn 's 6 ° B„=6T . 
(a) I Y A B „= 7T 



T 






I 




g'=2.5 


" Y>k 




x Xx - 

X 

x x 
x^ 


x x 


x^ 


^xx^ 1 


(b) 




v 6 n =8 T 



0.20 



E 0.15 

J5 

E 
0.10 



0.05 



FIG. 4: (color online) Experimentally deduced m*g* (open 
symbols), normalized to the free electron mass m e , as a func- 
tion of charge carrier concentration for SLG (a) and BLG (b). 
The individual data points were extracted from the total-field 
dependence of the SdH amplitudes corresponding to different 
Landau levels N — 2, ..., 10 and represent measurement for a 
constant magnetic field B n — 5 T, 6 T and 7 T for SLG and 
B n = 8 T for BLG. The error bars represent the standard 
least squares fitting errors, taking into account error bars of 
A (Fig. Esb. The dashed lines in (a) represent the calculated 
behavior of m* g* for different values of g* taking into account 
a 10% experimental uncertainty (shadowed areas). The blue 
crosses in (b) compare our data to the experimental cyclotron 
mass for BLG [17] multiplied by g* — 2.5. The inset shows 
the effective g-factor, extracted from the product m*g* in the 
main panel and the known cyclotron mass m* in SLG, as a 
function of Landau level index N. 



not symmetric for holes and electrons. Such a behav- 
ior is caused by an asymmetry of m* resulting from 
an asymmetric band structure of biased BLG, which 
was already observed experimentally in transport exper- 
iments [T7], cyclotron resonance [TB] and activation-gap 
measurements [5]. Applying the experimental cyclotron 
mass from Ref. d7| (depicted as blue crosses in FigJ4J) al- 
lows us to estimate g* to be about 2.5 for both electrons 
and holes which is, within experimental accuracy, reason- 
ably consistent with the g-factor enhancement observed 
in SLG. 

The observed enhancement of the effective spin- 
splitting compared to its free-electron value can be ex- 
plained by electron-electron interaction [TH] yielding an 
interaction-enhanced splitting between two spin levels 
within one Landau level [20] [2T] : 



g*VBB to t = gUsBtot + E° x (in - n t ). 



(5) 
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Here g — 2 is a free-electron g-factor, E® x is an exchange 
parameter, and n^ and n^ are the relative occupations of 
the two spin states of a given Landau level. 

For Gaussian shaped Landau levels with broadening 
F > g*[iBBtot, he. where the spin splitting is not yet re- 
solved, this relative occupation difference can be approx- 
imated using the Taylor expansion of the Gauss error 
function erf(g*/isi? tot /r): 
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is of the order of Coulomb interaction, i}g X ex 
y/R^ [H], and T oc y/B^ [25]. Therefore, the ratio E° x /T 
remains constant and the g-factor enhancement is indeed 
predicted to be constant as we observe experimentally. 
Using the experimentally found g* = 2.7 in Eq. ^ yields 
E° x = 130 K at 10 T when assuming r = 200 K g] [5] . 
For a completely spin polarized system, i.e. nj. — n^ = 1, 
one might then speculate that the exchange enhancement 
in the Eq. (pi) would be an order of magnitude larger than 
a single particle Zeeman energy at this particular field. 

Finally, we note, that the experimentally found en- 
hanced values of g* in graphene are close to those ob- 
served in transport experiments in graphite |23j . This 
may suggest that an exchange induced enhancement of 
g* is quite common for graphitic materials. In con- 
trast, no interaction-induced g-factor enhancement is ob- 
served using electron-spin resonance in graphene [24) and 
graphite [55] since these measurements are not sensitive 
to many body corrections [55]. Interestingly, measuring 
the Zeeman splitting of single-electron states in quantum 
dots, where no exchange enhancement of the (/-factor is 
expected, also yields g « 2 [27], albeit with a considerable 
experimental uncertainty. 

To conclude, we have experimentally measured and an- 
alyzed spin-splitting in SLG and BLG. We have shown 
that the product between the cyclotron mass m* and the 
effective g-factor g* increases with charge carrier concen- 
tration, as expected for a linear dispersion in SLG and a 
hyperbolic dispersion in BLG. Using the known concen- 
tration dependence of m* we found that g* in graphene 
is enhanced compared to the free-electron value and we 
attribute this to electron-electron interaction effects. 
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